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ABSTRACT ne 


The design and operation of a special purpose digital signal 
processing system is explained. The system consists of a Sienal 
analyzer, an interface betwsen the analyzer end the M.I.%, Conpat@= 
ible wime Sharing System, and appropriate commuter programs. With 
inputs from an instrweented shaft in the H.I.T. variable-pressure 


propeller tunnel, the system permits rapid spectral analysis of 


propeller generated forces and moments in six degrees of freedon. 
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1. INTRODUCTION 


1.1 BaCKGROUND. On July 1; 1968, the Department of Teva Archi- 
tecture and Marine Enginecring at M.I.T. submitted a proposal to the 
National Science Foundation to build a system that would permit une 
steady hydrodynamic measuroments in the variable-~pressure water tunnel 
operated by the department (1). nesearch ae means of reducing vib- 
ratory excitation to the hull of a ship depends upon the availability 
of methods to identify and measure periodic signals of known frequency, 
periodic signals of unknown frequency, and random signals. Of parti- 
cular interest are the six components of periedic forces and moments 
on the propeller shaft due to non-uniform flow about the propeller and 


the fact that the propeller has a finite number of blades. 


Ne 


The design and construction of a specially instrumented shaft 
which will enable simultaneous measurement of the instantaneous values 
of the three forces and three moments due to the propeller interaction 
is covered in the Naval Engineer's thesis by Charles 0. Horton (2). 
A matrix of miniature strain gauges and integrated circuit operational 
amplifiers built inside the shaft provides six signals in the 100 nv. 
to 1 volt range at an impedance of a few ohms. All the signals of 
interest will be below 1 khz. S&ach of the six forces and moments to 
be studied is a linear combination of these six signals, and the mate 
rix of coefficients for this computation is to be determined by dynae 
mic calibration of the shaft. A separate circuit in the shaft housing 


delivers a sharp synchronizing pulse once per revolution. Power to 
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the amplifiers and output signals from the amplifiers are transferred 


between the shaft and the outside via low noise slip rings. 


The design of a Simple digital data-processing system to rapidly 
analyze the signals received from the instrumented shaft is the sub- 


ject of this report. 


1.2 CHOICH OF SYSTEM. The principal requirement for the data. 
processing system were wide flexibility in types of analysis available 


to the exocrimenter, short turn-around time, and cost below $15,000. 


Originally it was felt that all signaleprocessing equipment 
should be self-contained and located in the propeller tunnel experi- 
mental area. Complete signal~processing systems which could correlate, 
Sienal average, end spectrum analyze turned out to be too expensive. 
Considered next was a method of recording data from the experinent, 
converting it to digital form, and then punching it on paper tape or 
cards for submittal to the M.1I.T. Information Processing Center for 
"batch" processing. -The long turn-around time made this prospect dis- 
couraging, The final choice was to use a conmercial digital signal 
analyzer in the experimental area, and to make an interface between 
the analyzer and the IF 709): commuter used in the M.1.T..Compatible 
Time Sharing System (CTSS). The analyzer would collect data from the 
exveriment, make the analog-to-digital conversion, and either store 
the data from a single sweep for analysis by an autocorrelating or 
crosscorrelating program in the 709), or make hundreds of sweeps in 


averaging mode, leaving in its memory a digitized signal of enhanced 





signal-to-noise ratio for use in a digital Fourier analysis. 


the signal analyzer chosen was the northern Scientific Model 550, 


with optional specifications as follows: 


Memory: 102), addresses; 10° counts per address 
Number of inputs; 8 (With NS-305 eight input expander) 
Readout: Teletype odsl 33ASKTC (includes tape 


punch and reader) 


the NS=550 mst the requirements of this project in the following 


specific ways: 


1. Because it covld read its output into a tcletype, which in turn 
can be commected to C1SS with a simple electrical interface to a vData~ 


Phone, the computer interfacing problein would not be comolicated. 


2. The system, in the worst case, should resolve the third harmonic 
of a 9~bladed propeller operating at 2000 repem. with all 8 inouts 
Simultaneously. With time multiplexing of the eight inputs, this 
requires a sampling frequency of at least J.) khz., or «a>.saindling 
period of no greater than 70 micreseconds, The NS=-550 can sample at 


intervals of 50 microseconds, 


3. ‘the total price of the analyzer and peripheral gear was about 


$12,000. 


whe electrical interface between the Teletype - machine and the 
Data-Phone was dosigned and built by the author. It is electronically 


very simole and requires no attention or adjustment, 
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The system is corpleted by a conversational program written in 
Fortran IV which receives typed inout from the exoerimenter and the 
Signel analyzer and returns a Fourier analysis of each of the six 
forces and moments. For an analysis involving all 102) memory loca- 
tions in the analyzer, the turn-around time is about 20 minutes, most 
of which is analyzer read-out time. The program listing and documen= 


tation is given in Appendix A. 


Hig. lis ablock diagram of the completed signal analyzing 


system. 





2. MODIFIGATLONS TO EQUIPMENT 


2.1 INPUr MODIFICATIONS. The Ns-550 analyzer has an input impedance 
of at least LOOK and switch-selected AC or DC coupling. The torque 
and thrust outputs from the shaft will have d.c. levels considerably 
higher than trrice the a.c. signal components, and consequently the 
baseline adjustments on the analyzer will be insufficient to remove 
them. Unfortunately, the time constant of the a.c. coupling is .01 
second,;: which is too small to pass low-frequency signals from the ex~ 
periment without significant attenuation and phase shift. This con- 
dition can be corrected by replacing the .1 pf. coupling capacitors in 
the eight input circuits with 10 pf. capacitors of equal voltage ra- 


tings. 


2.2 SWEEr MODIFICATIONS. With 8 inputs anda memory cycle time of 
50 microseconds, it takes .l) msec, to add one data point to each of 
the eight memory subdivisions. For simplicity, assume that each chan-e 
nel uses 125 memory locations. We desire to store data representing 
one complete propeller shaft revolution in each subdivision. The 
maximum period of one revolution is then 125 X .l) msec., which corres= 
Rens. to aminimum shaft speed of 1200 r.p.n. econ the speed 
does nothing more than reduce the number of data points corresponding 


to one period. 


Thus, for speeds of 1200 r.p.m. and a bone the time base of the 
NS-550 must be set aa 4 ms./point. Similarly, for spseds between 
600 and 1200 r.pen.; a setting of .6 ms./point is required. However, 
as supplied by the manufacturer, the NS-550 skivs over this range; 


LO 





going from .2 ms./point to 2 ms./point. 


The required modifications to the time base circuitry are not 
difficult, and circuit diagrams showing these changes can be obtained 


from Northern Scientific, © 


2.3 LINE FEED MODIFICATION>. Whenever CTSS is called up from a Tele~ 
type terminal, the computer begins by signaling the automatic answer- 
back device on the Teletype. The codes which are returned identify 
the particular Teletype machine and tell the computer if it mst sun- 
ply separate line feed signals whenever a carriage return is either 
sent or received. The two codes are identical except that in one a 
letter "C" appears in the two places which contain the letter "D" in 
the other code. If the answer-back device is set up for the "D" codes; 
the computer supplies a separate line feed. If it is set up for the 
“C" code, the commuter does not supply the separate line feed, but age 
sumes that whenever a carriage return is either sent or received, the 


Teletype machine will line feed on its own. 


The "D" code works satisfactorily for ajl input and output exe 
cept for signal analyzer readout, because the signal analyzer produces 
a line feed of its own, and does not wait for the receipt of the come 
puter-generated line feed before continuing. This causes confusion 


and loss of data. 


Using the "C" code solves the above problem and introduces ano- 


ther: the operator must type his own carriage returns when typing 
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input to the computer, and must manually advance the paper for each 


Successive line of data returned by the computer. 


The solution to this problem is to set up the function box in the 
teletype machine for "Suppress Line Feed" and "New Line, Line Feed," 
in accordance with the technical manual supplied with the machine (3). 
the “Suppress Line Feed" feature ignores all line feed signals re- 
ceived, and the "New Line, Line reed" function automatically produ- 
ces a line feed whenever a carriage return Signal is either transmit- 


ted or received, vf course, the "C" answer--back code is used. 
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3. ELECTRICaL INTERFACE BETWEEN TELETYPE AND DaTa-PHONE 


3.1. ‘TELETY+e SIGNALS. she 33ASRTC Teletype that Northorn Scientific 
supplics as an optional readout device for the NS-550 comes wired by 
them for three modes of operation: “Tape to Type," "Analyzer kead 
fn;" and "Analyzer Read Out." In the first mode, the tape reader sup- 
plies data to the printer. In the second mode, the tape rmder sup- 
plies data only to the memory of the analyzer. Tne last made has the 
printer (and tape punch, if desired) print (punch) out the data in the 
analyzer's memory. . In all these modes, the data is transferred in 
parallel to and from the analyzer. The parallel-to-serial conversion 
of the code for each symbol (required for transmission via the Data- 


Phone) is accomplished electromechanically by the printer mechanism. 


The Teletype itself has provision for both "local" and "line"! 
operation. In local mode, an internal power supply provides the 60 
mae loop current required to operate the selector magnet system. When-~ 
ever the machine is being used for purposes not requiring a connection 
to the conputer, it should be operated in local mode. In line mode, 
an external power supply is required, and it is connected to terminals 
6 and 7 of terminal strip X. (For details of the ‘feletype wiring, see 
reference (1). If the current from the external supply is modulated 
on and off by the ASCII code for a letter, the machine will print the 
letter (and punch it as well, if desired). Alternatively, if a key on 
the keyboard is deppressed, the Teletype itself modulates the line 


current with the ASCII code. This also happens when a symbol is 
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printed during operation in “Tape to Type" or in “Analyzer Read Out" 


mode e 


3.2 DATA-rHON SIGNALS. The Data-Phone transmits and receives all 
data serially, in the form of pulses of (nominally) plus and minus 15 
volts. Full details are given in the manufacturer's interface speci- 


fication (5). 


3.3 THE INTERPACE. The interface translates the modulated line cure 
rent from the Teletype to positive and negative voltages for the Data- 
Phone, and vice Verse Fig. 2 is a schematic diagram of the unit. 
The interface consists of a 60 ma. loop supply, at10 volt supply, and 
two transistorized switches. The incoming signal appears at terminal 
3 of the "customer's equioment" connector. Whenever the "Mark" (nega- 
tive) sienal is received, the loop current is switched on by Q1. ~The 
liSpace" (positive) signal causes Ql to switch off. The germanium 
diode connected to the base of Q2 serves two purposes. First, its 
back resistance is sufficient to assure that Q2 is cut off whenever 
the Data-Phone is disconnected, Second, it prevents damage to the 
base-to-emitter junction in the event the Date—-Phone produces a nega- 
tive voltage greater than -10 ah Gensel 3. (It may go as high as 


-25, according to the interface specification.) 


Q3 and Qi supply the oubgoing signal, The voltage divided con- 
nected to the bases of Q3 and Q; is chosen so that when the Teletype 


is marking (terminals 6 and 7 of strip X shorted) a negative voltage 


a 
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at low impedance is presenved at pin 2, and when the Teletype is spa- 
cing, a corresponding positive voltage at-low impedance appears-at*pin 


26 


The continuovs positive voltage connected to pin 20 enables the 


Data-Phone to enter and remain in the "Data! mode, 
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lh. PROGRAMMING REQUIREMENTS 


h.1 OUTPUT FORMAT FOR THE NS-550. <A typical line of output from the 
NS-550 consists of a line feed (LF), a carriage return (CR), and eight 
6-digit numbers, each separated by a single space from the preceeding. 
If each memory location contains only 5 decades (which is the ease On 
the NS=-550 used in these exneriments), the first digit of each number 
is always zero. The analyzer never prints e plus or minus sign. With 
the modification described in section 2.3 of this report, the Teletype 
' ignores the line feed, but automatically produces one when it receives 


the carriage return. 


The initial CR and LF is interpreted by the computer as a blank 
card image, and in order to avoid having the first 8 locations of the 
input array set equal to zero, the program must be written to read the 


blank image before reading the first 8 data points. 


The following program segment will successfully read JMAX points 


from the signal analyzer: 


READ 100 
100 FORMAT (1X) 

DO 10 N=l,JMAxX, 8 

NA= N¥7 
10 READ 101, (INDATA(NZ) ,NZ*N,NA) 
101 FORMAT (8(16,1X)) 
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4.2 AuLOCATION OF STORAGE IN THE NS-550. ‘he 102) memory locations 
in the analyzer are divided equally ‘among the input channels being 
used. ‘The first and last locations in each block of memory are always 
set to zero, except for the first Jocation in the first block, which 
always contains the number of sweeps made in the experiment. vata is 
always read out sequentially by memory address number. Data from Chane 
ne] A appears in the first block of. memory, data from vhannel B ap- 
pears in the second block, and so on. A computed GO TO statement with 
its argument equal to the number of input channels can be used to as- 
Sign the data received from the analyzer readout to arrays indexed by 


channel number (See Appendix A). 


4.3 INPU?S sUrPLIED BY THE EXPERIMENIER. Because time-sharing allows 

conversational interaction between the program and the experimenter, 

the program can ask for the values of various constants and receive 

the typed responses directly. The following sets of constants are re- 

Giaared : 

1. The number of input channels used in the experiment 

2. The settings of the input attenuators on the analyzer 

3. The matrix of coefficients relating the actual forces and moments 

to the Signals appearing at the input of the SHEESH 

lh. If the experiment involved signal averaging, and the readout did 
not begin with the first address of memory: the number of sweeps 


made by the analyzer 


5. if the program does a Fourier analysis: the number of data pdints 
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in a channel which represent one conplete shaft revolution, and 
the number of harmonics to be computed 

If the program computes a correlation funetions the range of the 
delay time (the argument of the correlation function) 

If the program crosscorrelates: the values of the second function 


or* an algorithm for computing these values 


the Fourier analyzing program in Appendix A has numerous examples 


of conversational inout-output routines. 


ay 





5. OPERATING INSTRUCTIONS FOR A SIX-DIMENSIONAL ANALYoIS 


It is assumed in these instructions that the reader is familiar 
with Mr. Horton's thesis (2) which describes the overall design of the 
complete experimental system, the operation of the NS-550, which is 
covered in the instruction manual (6), and the UTSS procedures, which 


are covered in the CTSS User's Manual (7). 
PROCEDURE 


5.1 Push the STOP button on the NS~550. 
5.2 Set: DECADE SELECT to 10” 
READ MODE to CRI 
ADD-SUB to ADD 
MEMORY GROUP to 1/1 
VIME BASE to 2 ms./point 

NO. OF SWEEPS TOo° 

DELAY TIME to lms. 

PRE DELAY,’ POST DELAY to OFF 

ADDRESS SELECT to OFF 

VERTICAL OFFSET to OFF 

CONVERSION GAIN tu saVE 

SWEEP MODE to AVE 

"RIGGER to RECURRENT 

VIEW INPUT<-v EW MEMORY to VIEW MeMURY 


INPUT SELECT to nxt 
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5 ly 


5.5 


Sor 


5.8 


Push START READOUT. «a pattern of some type should appear on the 
CRT. if nothing appears , carefully adjust the INTENSITY, HORI- 
ZONTAL, and VERTICAL controls until something does, Always keep 
the INTENSITY control set to the minimum level which produces a 
visible trace or spot on the CRY. This is especially imoortant 
when only a single spot appears on the CRT. 

Push the ERAS# button. The CRT pattern should immediately degen- 
erate to a straight line. Adjust the VERTICAL POSITION control 
so this Jine appears at the very botton SageerOl Lhe sGiime @Ce. 
Adjust the HORIZONTAL POSIT:«ON and HORIZUNTAL EXPAND controls so 
that this line extends slightly less than the width of the CRT 
face. 

Push STOP, Then push START MEASURE. Let the unit run for about 
10 seconds, then push STOP, followed by START READOUT. This cas~— 
ses a small amount of noise to be read into the memory. 

Set the DECADE SELECT to ilo and anstercne that 10 horizontal lines 
appear on the CRT. Using the VERTICAL EXPAND and VERTICAL POSI- 
TION controls, make these 10 lines coincide with 10 of the 11 
scribed reticle lines. This completes calibration of the CRT. 
Push STOP, Set DECADE SELECT to 10°. Set the INPUL ATTENUATURS 
for all l, inputs of the NS-550 ‘and all i inputs of the NS-305 
to l volt. Set all 8 of the AC-DC switches to AC. 

Connect the six strain-gauge outputs from the shaft to channels 
A through F. Connect the synchronizing pulse cable from the 


shaft to the positive input of Ghannel H and also to the EXT 
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5.10 


5.11 


Sale 


5.13 


5.Uh 


TRIGGER connector. 

Start up the propeller and set the shaft r.p.m. to the value that 
will be used in the emsriment. Set the TIME BASE to .) ms./ 
point if the r.pem. is above 1200. Set the VIE BASE to .8 ms./ 
point if the speed is between 600 and 1200 rep.m. Set the TIME 
BASE of the NS-305 to the same value that is set on the NS-550. 
Set Vibw INPUT to A. Set TRIGGER to +EXT, and push START MEASURE. 
Adjust the TRIGGER SENSITIVITY for stable triggering of the hori- 
zontal sweep, as indicated by a stéady pattern on the CRT that 
does not move left or right. Center the pattern vertically by 
using the Channel A BASnLING ADJUST. If the pattern goes off the 
top or bottom of the CRT, set the Channel A INtUT ATTENUATOR to 
10 volts. If the pattern is very small, set the Channel A INPUT 
ATTENUATUR to 0.1 volt. mt 
Switch VW INvUL GO 7b. Using the Channel B BASELINS ADJUST and 
INPUT ATTENUATOR, position and vertically scale the pattern as 
before. nepeat for Channels C and D. 

Sev the DEC.DE S#LeCl’ to LOG, and Vii INPUT~VinW MEMORY to VIEW 
MeMOnY (shis permits viewing of the four inputs of the NS~305). 
Set VIEW INPUT on the No=-305 to BE. using the Gheimel EB INPUT AT} 
JENUATOR and BASELINE ADJUST, position and scale the pattern as 
before, Kepeat for Chammel F. 

Set the ViseW IN. Uru to G. Genter the line using the Uhannel G 
BASELING ADJUST. 


Set VIEW INPUT to H. whe pattern should be a straight horizontal 
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5.15 


5.16 


5.17 


S adls 


Soule 


line except for a pulse near the right-hand side. With the BASKH- 
LINE adjust for Channel H, position the straight line about 19% 
up from the bottom, With the INPUs ATYENUATOR, make the pulse 
heignt at leest two divisions on the reticle. 

Push STOP, Set DECADE SELECT to 10, Push STANT READOUT. Push 
KRASE, sush STOP. This completes the adjustments, 

To begin the measurement, push START MBASURE, Initially the CRT 
will show only a horizontal line ab Zoro. As the experiment con- 
tinves, vertical deflections will appear and grow. Set DECADE 
SELECT to 10°. ‘Continue the experiment until the position of any 
portion of the trace approaches the top of the CRT. ‘hen push 
STOP. (Alternative: Set the nO. OF SWEEPS control to 800 before 
pushing START MEASURE, ‘the analyzer will make 600 sweeps and 
stop automatically.) ‘this completes the measurement phase of the | 
experiment . 

Turn on the electrical interface between the Teletype and the 
Data=Phone., Turn on the Teletype in "Line mode, Set the funce~ 
tion switch to "ANALYZER READ OUT." Dial up the computer and 
LOGIN. Then type LOADGO PROP, The program will begin by askinz 
ce a number of constants. ‘fhe following steps give the details. 
"TYPE NUMBER OF INrUT CHANNELS," Format is Il. Type in the num- 
ber six. (The sync input on Channel H is not. counted as an ine 
put.) 

NTVPi INPUT ATTENUsTOn SETTINGS." Format is 6F7.2 Type in the 


settings of the six INPUT ATTENUATORS, A through F. 
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5.20 "TY.E NUMBER OF OUTPUT CHANNWLS DESIRED." Format is Il, Type 6, 


5.21 "TYPE NUMBsR OF HARMONICS DESIRED." Format is I2. The basic 


562 


5.2 


5.2 


& 
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MN 


period of analysis is one shaft revolution, EXAMrLE: To obtain 
analysis up through the third harmonic of the blade rate for a 5- 
bladed propeller, 15 harmonics are required. The maximum number 
that can be typed in is 50. 

"ENTER THE N-OUT BY N-IN CALIBRATION MATRIX." The response must 
be the 6 X 6 strain-gauge calibration matrix C,,s Where k, the 
row subscript, is the output channel number, andi, the column 
subscript, is aie input channel number, The dimensions of each 
entry must be either ft.-lbf./volt or 1lbf./volt,; and the computer 
will read six successive rows of six numbers each on a 61'7.2 fore 
mat. 

Wty NUMBER OF DATA POINYS TO BS READ." Format is Ih. Type in 
the number 102), 

"HERE IS HOW I HAVE YOUR SETTINGS-" The computer will read out 
all the values as it has them, and gives self-explanatory in- 


structions on wnat to do if there is an error anywhere. 


NOYE: In the next two instructions, the program mentions switch- 
ing to either manual or automatic lines feed. If the modifica- 
tions described in section 2.3 of this report have been made to 
the Teletype, the line feed instructions should be deleted from 


the program. 


5.25 "I AM READY TO ACCEPT DATA FROM THE SIGNAL AVERAGER." If a papor 
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5.27 


5.20 


tape ‘of the readout is desired, turn on the tape punch and hit 
the space bar on the keyboard about 10 times to make a leader on 
the tape, but do not type a carriage return or a line feed. Push 
STOP on the NS-550., Set the READ MODE to TYPE and push START 
READOUT. The readout will take about 15 minutes. 

USWITCH TO AUTOMATIC LINE FEED AND TYPE CARRIAGs, RETURN." See 
the note above. If this instruction and the read Satoneneens 
mediately following it in the prograi (Appendix A) have not been 
deleted, simoly hit the carriage return to continue. 

“TYPE NUMBER OF DATA POINTS GOMPRISING OWE CUMPLErE PERIOD." In 
the event the program was not able to detect the sync pulse in 
Channel 8, it will make the above request, If this happens, look 
at the last 16 rows of the printed readout. Observe that the 
sixteenth row from the bottom begins with 000000. Number these 
rows 1 through 16 starting with the one beginning with zeroes. 
Then find the discontinuity in the printed data values. If R is 
the row and C is the column where the discontinuity appears, the 
number to types is 8(R-1L)+(C-1). The format is Th. 

UREaD OUL NUMBER OF paSSES." If the comouter receives a zero for 
re first number read in rom he analyzer, it will make the 
above request. Push STOP on ae Ns=-550. Set MeMoRY GROUP to 1/1 
and turn the ADDRESS SELECT switch OF, Set READ MODE to ‘CRT, 
Push START READOUT and then push STUP, Set READ MODE to TYPE and 
push START READOUT. After the first number is printed and before 
the eighth is printed, push S!OP, set READ MODE to CRT, and then 


hit the carriage return. 


N) 
AR 


{ 
5.29 The computer will now read out by channel number and harmonic 


5.30 


number the amplitude in lbf. or ft.-lbf and the phase in dege 
ree3 of each harmonic, The zeroth harmonic is the average, and 
it has no significance in this typos of expcriment, since the in- 
puts are all a.c. cowoled. ‘the channel numbers can be related to 
a specific force or moment by a brief consideration of the compus 
tation innplied by the 6 x 6 calibration matrix. 

"THIS COMPLETES THE ANALYSIS. TYPE 1 10 START AGAIN, O vO STOP." 


Tf a "1" is typed, the program will return to step 5.25. 
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APPENDIX As "PROP," A COMPUTER PROGRAM FOR FOUKRTER ANALYSIS 


The program described 3n this section embodies the principles 
outlined in section l} and is the progrem used for the 6-dimansional 
analysis of section 5. The listing begins on page 32. PROP can be 


divided into 5 parts: 


1. TYPED INPUTS: PROP requests values from the experimenter and pro- 
vides an opportunity to start again if errors have occurred, (Cards 


1-37) 


2. ANALYZER READOUT: The program reads the analyzer output into the 
array INDATA, and then assigns the values to appropriate blocks in the 
array INPUT, according to the rule stated in section 4.2. If 5 or 
more snpwt channels are being used, the program scans the contents of 
Channel 8, looking for a discontinuity of at least 20,000 counts. If 
it finds such a discontinuity, it sets the number of data points rep~ 
resenting one comolete shaft revolution equal to the location of the 
discontinuity. If the discontinuity is not found, or if the number of 
inputs is } or less, PROP requests the experimenter to type in the 


correct value. (Cards 38-87) 


3. MATRIX MULTIPLICATION: The program converts strain-gauze signals 


into forces and moments. (Cards 88-98) 


ht. FOURTER ANALYSIS: PROP continues by calculating the coefficients 
for the Fourier sine and cosins series for each force and moment. 


(Cards 99-116) 
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5. OUTrvUr: The program determines the number of signal analyzer pas- 
ses (sweeps), calculates the magnitude and phase of each harmonic for 
each output charmel, prints the results, and provides the experimenter 


the opportunity to read in new data for the next analysis. (Cards 117- 


152) 


The following 


is a chart of the variables used in the progran. 


VARIABLE NAME RANGE Ov VALUES DESCRIN TION 
A(K5L) 1bf./100 counts or Fouricr sine coefficient for 
ft,-1bf./100 counts the Lth harmonic in the Kth 
output channel 
AMP fG,—Lbr. or bf. The magnitude of a harmonic 
AVE” ft.-lbf. or 1bf. whe average value of the 
Signal in an output channel 
Bik, L) 1bf./100 counts or Fourier cosine coofficient 
; ft.-lbf./100 counts for the Lth harmonic in the 
Kth output channel 
C 2T/RDPC Coefficient used in Fourier 
series calculations 
CALIB(K,<) tpi volt or The k,i entry in the strain- 
ion, / TOIL gauge calibration matrix 
FULSCA(T) 0025-100, volts per Input attemator setting of 
100 counts input channel i 
i 1-8 Index giving input channel 
number 
IMAX 168 Number of input channels be- 
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INDATA(N) 


InPUL(I,J) 


Jy JP, JPO, @UC o 


JMAX 
K 


L 


M 


NA 


NDPC 


NOUTCH 

NX 3 NY 

NZ 
OUTPUT(K,J) 


PASSES 


PHASE 


0-99999 
0=99999 
1L--1.022 


1-102) 
1-8 


1-50 


THe) 


6, 16, 2h, etc, 


LelO22 


1-8 


Qe] 
L102) 


1b./100 counts or 
ft.-lbf./100 counts 


Le-99772 


90, to +2/0. desrees 


390 


Yhe Nth data point read in 


from the signal analyzer 


The Jth data point in input 


channel. i 


Index giving data point mime- 
ber 


Number of data points read 


Ovtpyt channel index number 


Index of harmonic being cale 
culated 


Number of harmonics the pro~ 


gram is to calculate for 


each outout channel 
Index vwsed in reading INDATA 


Number of data points per 
channe] which represent one 
conplete period of the inout 


Signal 
Number of output channels 


Branching variables 
Index used in reading INDATA 


The Jth data point in output 
channel k 
The number of sweeps made by 


the 


experiment 


Sienal analyzer in an 


The phase of. a harmonic in 
reference to a sine function 


of the same frequency 





RIN 0.-99999, Real value of 1LNPUT(I,J) 


RJ ager Real value of J 
scos(J)"" 7 “1.0 to 1.0 Stored value of a cosine 
SSIN(J)” ia ye) we) 150 Stored value of a sine 


ee a ee a Oe ee rE a 8 Ge eh ee OF OS ee 


“AVE is printed as the Oth harmonic in the output format. 


3% : 
“PHASE is initially calewated in radiens from the arctengent lib- 


ree 
we 


rary routine, and then is converted to degrees, taking account of 
the quadrant, 


“SSIN and SCOS are used to avoid repeated computations of the sine 
and cosine values needed for the Fourier analysis of a given har-: 


monic in diiierent output channels. 
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APPENDIX B: NOTES 


1. Representing the low-impedance source as an ideal voltage source, 
the transfer function for the capacitor counled inovt is 


JzmRcFf 


H (f) = (E7ER RCE 


A straightforward calewlation shows that the phase of H(f) is 
op (7) 
?=tan Zit RCF 


and the magnitude is 


[ 


HAs Vie aaa 


fee poe shaft r.p.m. set to 950, f is 16 hay. and with RC .01 sec., 
& Cpe 
anRcfad., $2 45° and [HOLS F= 


Increasing the value of RC to 1 second permits reducing the r.p.m. to 
only 96. Under these conditions, the phase shift is only 5.7 degrees 


and the magnitude is down by a factor of only .005. 


2. Northern Scientific, Inc., 2551 West Beltline, P.O. Box 66, Mid- 


dleton, Wisconsin, 53562, Phone (608) 835-6511, 


3. In the computer program of Appendix A; each entry in the calibra= 
tion matrix is assumed to be independent of frequency. If dynamic 
calibration of the shaft shows that this asswmtion is not valid, the 
program will have to be rewritten to do a Fourier analysis of the in-~ 
coming data before doing the conversion of strain-gauge signals to 
force and moment cutputs. wzhen each harmonic will have to be converte 
ed separately in the program, vsing values for the calibration matrix 


entries that are appropriate for the frequency of the harmonic, 
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